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SECTION  I 
INTRODUCTION 


The  problem  of  electromagnetic  scattering  by  chaff  clouds  can  be 
treated  in  two  distinct  fashions,  namely,  deterministic  and  nondeter 
ministic  or  statistical.  As  might  be  expected,  each  has  its  own  merits 
and  drawbacks.  Deterministic  methods  predict  the  exact  outcome  in  any 
given  situation  and  are  generally  quite  involved.  In  the  case  of  chaff 
clouds,  there  are  two  separate  parts  to  the  problem,  the  aerodynamic 
and  the  electromagnetic.  First  one  computes  the  orientations  and  posi- 
tions of  the  individual  dipoles  in  the  chaff  cloud  using  a generalized 
six-degree-of-freedos-equations  program  and  then  the  scattering  be- 
havior due  to  plane  wave  excitation  determined  by  means  of  well- 
documented  standard  techniques.  In  each  case,  a digital  computer  is 
essential;  but  for  the  digital  computer,  such  calculations  would  be 
unthinkable.  Furthermore,  the  effects  of  environment  such  as  wind 
shear  and  turbulence  can  be  included  as  excitation  parameters  in  the 
aerodynamics  calculations.  Deterministic  methods  in  the  present  case 
are  time-consuming  and  very  expensive,  and  in  view  of  the  basic  limita- 
tions of  computers  such  as  finite  memory,  etc.,  the  number  of  dipoles 
that  can  be  handled  cannot  possibly  exceed  a few  hundred.  These  tech- 
niques are  being  investigated  by  the  Air  Force  Avionics  Laboratory 
through  other  contractual  efforts  and  will  be  reported  later. 

Statistical  methods,  on  the  other  hand,  are  not  concerned  with  any 
one  particular  situation}  they  predict  in  a probabilistic  fashion  what 
might  happen  under  a given  set  of  circumstances.  The  methods  are  quite 
general  and  apply  to  a variety  of  problems  occurring  in  physics  and 
engineering.  While  deterministic  methods  are  severely  limited  by  the 
number  of  dipoles  that  can  be  handled  in  the  caec  of  chaff,  there  ie 
no  such  restriction  with  statistical  methods}  in  fact,  the  larger  the 
number  of  dipoles,  the  more  accurate  the  predictions  become.  This  re- 
port treats  the  problem  o.  electromagnetic  scattering  from  chaff  clouds 
in  a systematic  and  quar  tit&tivc  manner  using  statistical  methods.  The 
eater ial  presented  tsere  has  been  gathered  from  different  sources  and 
put  into  self-contained  form.  Certain  mathematical  derivations  have 
bean  simplified  considerably. 
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The  plan  of  this  report  is  as  follows.  First,  basic  material  on 
probability  and  random  processes  is  introduced.  Then,  first  and  second 
probability  densities  of  chaff  cloud  scattering  are  derived  starting 
from  first  principles.  Relevant  averages  and  auto-correlation  are 
obtained.  The  physical  significance  of  each  random  function  is  ex- 
plained with  illustrations.  A self-consistent  mathematical  model  for 
chaff  cloud  scattering  is  developed.  This  will  be  a basis  for  further 
analytical  studies  in  such  important  areas  as  the  effects  of  chaff  echo 
fluctuations  on  continuous  wave,  pulse  doppler,  and  MTI  radars. 
Furthermore,  we  hope  to  combine  statistical  and  deterministic  methods 
in  a judicious  manner  so  as  to  be  able  to  predict  chaff  cloud  behavior 
more  accurately. 


SECTION  i: 
PROBABILITY  THEORY 


2.1  DEFINITIONS 

A random  variable  (r.v. ) that  is  a function  of  time  is  called  a 
random  process.  Let  us  denote  the  r.v.  by  x(t)  and  che  value  attained 
at  an  instant  of  time  t^  by 

xk“x(*k)  (2-v 

Because  of  the  random  nature,  it  is  meaningless  to  talk  about  the 
value  attained  at  a particular  instant  of  time  or  the  values  observed 
over  a period  of  time.  Then  how  does  one  handle  the  problem?  The 
answer  is,  of  course,  by  using  the  notion  of  probability  which  lends 
to  the  precise  definitions  of  certain  distributions  and  averages  which 
can  be  predicted  and  observed  with  some  measure  of  confidence.  The 
first  and  second  order  probability  density  functions  (PDF)  are  defined 
by 

P1  (x;t)  = Pr°bability  of  finding  x between 

x and  x*dx  at  time  t 

P2  ( xi *x2 • ^ 1 * 1 2 ) dXjdx^  “ J°int  probability  of 

finding  a pair  of  values 
x in  the  ranges  {x^,x+dXj) 

at  time  t ^ and  (x?,x2*dx,) 

at  time  t2 


These  definitions  can  be  extended  to  still  higher  order  PDF,  but  they 
will  not  be  needed  in  our  investigations.  For  convenience  of  writing, 
the  differential  g1m*o  .s  will  not  be  carried  along  with  p,  and  p_, 
but  should  always  be  understood.  In  the  definition  of  p„,  it  should 
be  understood  that  both  x^  and  x^  are  random  variables  with  the  same 
distribution  p^.  These  two  variables  ore  considered  statistically 
independent  or  simply  independent  if 

(xl,xl 'll,tl)  *>l(x2*t2)  (2-2) 
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- P2  ( xi  *ti  tx2*t:2  ) ?!  ( x2 ' t2  ) 


and  and  p2  fulfill  the  relations 


/ Pi  ( Vh)  dxi 


ss  1 


(2.4) 


I ?2  (3t1.*,it1.t2)  dx2  - Pl  (x1,t1  ) (2.5) 

f P2  (x2;t2lxlill)  Ux2  " 1 


In  relation  to  the  independent  variable  time  t,  random  processes 
are  divided  into  two  categories.  These  are  stationary  and  nonstationary 
processes  according  as  to  whether  the  statistics  are  independent  of  or 
dependent  upon  t.  The  latter  processes  are  extremely  complex  and  will 
not  be  considered  here.  Stationary  processes  are  further  subdivided 
into  strictly  stationary  and  wide- sense  stationary.  For  our  purposes , 
it  suffices  to  consider  only  wide-sense  stationary  processes.  What  this 
means  is  the  following.  The  first  order  PDF  is  indop  ndont  of  time  and 

the  second  order  POF  depends  only  c-r  the  difference  t * T.  Hence 

*»  1 
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Such,  of  course,  is  not  always  the  case.  To  handle  the  dependent  case, 
the  notion  of  conditional  probability  is  introduced.  This  is  denoted 
by  P2*x2't2lXl'V  which  9^v®*  the  probability  of  finding  x2  in  the 
range  (x2,x2  + dx2)  at  time  t2  given  that  x - x^  at  time  t^.  A vertical 
bar  separates  the  two  sets  of  variables  and  the  variables  appearing  on 
the  right  side  have  already  occurred  and  are  considered  no  longer 
random.  One  has,  by  definition 

p2  (xi»x2;tl*t2)  * p2  (x2,t2lxl»tl)  P1  (xl'fcl) 

(2.3) 

“ p2  ( X1  »tl  lx2,t:2  ) P1  ( x2  ft2  ) 


and  p^  and  p2  fulfill  the  relations 


/ i’i  ( xi;ti)  dxi  " 1 (2-4) 

/ p2  (xl,x2,tl’t2)  dx2  * P1  (xl’tl)  C2-5) 

/ P2  (x2:t2lXl!tl)  dx2  " 1 (2,6) 


In  relation  to  the  independent  variable  time  t,  random  processes 
are  divided  into  two  categories.  These  are  stationary  and  nonstationary 
processes  according  as  to  whether  the  statistics  are  independent  of  or 
dependent  upon  t.  The  latter  processes  are  extremely  complex  and  will 
not  be  considered  here.  Stationary  processes  are  further  subdivided 
into  strictly  stationary  and  wide>sense  stationary.  For  our  purposes, 
it  suffices  to  consider  only  wide-sense  stationary  processes.  What  this 
means  is  the  following.  The  first  order  PDF  is  independent  of  time  and 
the  second  order  PDF  depends  only  on  the  difference  tj-t^  « x.  Hence 
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(2.7) 


Pi  (x;t)  - Pl  (x) 
p2  (xi»x2>tl>t2)  = ^2  (xl,x2;t)  (2.8) 

for  wide-sense  stationary  processes. 

The  notion  of  homegeneity  in  time  is  sometimes  employed  to  de- 
scribe the  foregoing  random  process.  One  should,  of  course,  not  lose 
sight  of  the  fact  that  x is  still  a function  of  time.  From  now  on, 
subscripts  denoting  the  order  of  the  PDF  will  also  be  omitted  for  con- 
venience of  writing. 

The  cumulative  probability  denoted  by  P(x)  which  gives  the  proba- 
bility that  < x < a is  defined  by 

P(a)  = Prob(-“<x<a)  = /a  p(x)dx  (2.9) 

- 00 

from  which  it  follows  that 

p(x)  3 ^Lp(x)  (2.10) 

The  complementary  function 

Pc(x)  3 1 - P(x)  (2.11) 

is  also  commonly  used. 

The  expectation  or  expected  value  of  a function  f (x)  is  defined  as 

E(f)  = / f(x)p(x)dx  (2.12) 

with  the  most  important  ones  being 

E(x)  3 / xp(x)dx  (2.13) 

E (x2  ) 3 / x2p(x)dx 


(2.14) 


The  variance  and  standard  deviation  (SD)  of  x become 

Var  = D(x)  = E (x2)  - E2(x) 

SD  = o (x)  = /D(x) 

If  x represents  voltage  for  instance,  E(x)  is  the  D.C.  component, 

E(x2)  the  mean  square  and  D(x)  the  A.C.  component.  In  relation  to  the 
second  order  PDF  the  most  significant  quantity  of  interest  is  the  auto- 
correlation function  B(T)  defined  as 


(2.15)  | 

i 


B(t)  = E (xxx2)  = E(x(t)x(t+  c) ) 

= //  xLx2  p (x1,x2;t\  dxtdx2  (2.16) 


It  may  be  noted  that 

b(o)  - E (x2)  (2.17) 

and  for  large  time  lags,  x^  and  x^  will  be  uncorrelated  so  that 

B(~)  = E2(x)  (2.18) 

which  is  again  the  D.C.  component. 

The  auto-covariance  K(T ) and  the  normalized  version  R(t)  are  de- 
fined by 

K(t)  =*  B(t)-B(«)  (2.19) 

R(t)  - (2.20) 


It  may  be  noted  that  -l  < R(t)  < 1.  Since  the  correlation  function 
gives  the  correlation  between  x(t)  and  x(t  + t) , the  more  rapidly  x(t) 
changes  with  time,  the  more  rapidly  R(t)  decreases  from  its  maximum 
value  of  unity.  This  decrease  may  be  characterized  by  a correlation 
time  t defined  by 
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(lo) 


i 


i 

\ 


i 

% 


1/e 


(2.21) 


where  e is  the  base  of  natural  logarithim.  The  foregoing  gives  us  a 
clue  that  R(t)  and  the  frequency  spectrum  must  somehow  be  connected. 
This  indeed  is  the  case  and  given  by  the  Wiener-’. Chin tchine  theorem. 

a 

4>  O)  = Jo  R(tj  cos  bit  dt  (2.22) 

00 

R ( t ) » j f i (b) ) COS  tb)  db)  (2.23) 

It  is  noted  that  the  frequency  spectrum  $ rnd  R fort  a Fouiier  tr;  is- 
form  pair.  It  will  not  always  be  possible  to  integrate  (2.22),  and  in 
such  cases,  it  is  customary  to  take  b)  v - 2Ti/tQ. 

Now  let  us  introduce  the  concept  of  time  averages.  Suppose  there 
are  great  numbers  of  identical  radio  receivers  (ensemble)  turned  on 
simultaneously.  Let  us  also  assume  that  the  transients  have  died  down 
and  steady  conditions  have  been  reached.  The  noise  output  voltages  of 
all  the  receivers  are  recorded  over  a long  period  of  time  T.  At  a 
definite  time  t^,  we  take  the  voltages  x^(t^),  x^(tj),...,  compute 
the  average  and  the  probability  density  function.  This  average  is 
called  statistical  or  ensemble  average  and  this  is  what  we  have  con- 
sidered thus  far.  Stationarity  in  this  context  means  the  statistics 
are  the  same  regardless  of  the  value  of  t.  We  might  just  as  well  take 
the  output  of  a single  receiver  and  define  a time  average  in  the 
customary  manner  (denoted  by  a over  bar) 

STtT  - Jr  J.t  ’■(*><“  <2-24) 
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, From  an  experimental  veiwpoint,  it  is  much  easier  and  more  con- 

venient to  measure  time  averages.  Naturally  we  would  like  to  know  the 
relation  between  time  average  x(t)  and  ensemble  average  E(x).  Under  the 

• so-called  ergodic  hypothesis,  these  two  are  equal.  This  identity  will 

be  invoked  here  as  a basis  for  comparing  theory  and  experiment.  The 

! 

j time  average  auto-correlation  is  defined  by 
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B(t)  = x ( t ) x ( t + T j 
T 


T->»  IT  /-T  x^  x (t+^dt 


(2.25) 


taking. a single  realization  of  the  random  process. 


Under  the  ergodic  hypothesis  for  correlation  functions 

B(x)  - B(t) 

2.2  FUNCTIONS  OF  RANDOM  VARIABLES 


(2.26) 


Suppose 

y = f(x) 

and  we  want  to  determine  the  PDF  of  y.  If  the  inverse  function 


(2.27) 


x*  f 1 (y) 


(2.28) 


is  single  valued  <±en 


py(y) 


Px[f  1(y)]  ls£l  (2.29) 


layl 


where  subscripts  are  used  to  distinguish  different  functions.  In  case 
of  multiple  values,  we  first  define  single  valued  branches 


X1  “ *x2  = f2  W 


(2.30) 


and  get  the  more  general  formula 


Vy)  ‘ p*  (x0 


dxj! 


ar 


+Px  (x2  ) 


dx7 

ar 


(2.31) 


These  idear  are  easily  extended  to  functions  of  several  random 
v^riaoles.  For  example,  if  x and  y are  random  variables  and 
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u = u(x,y),  V = v(x,y) 


(2.32) 


then  in  terms  of  the  joint  density  of  x and  y 

Puv(u,v)  - Pxy[x(u,v)  ,y(u,v)]  (2.33) 

where  we  have  assumed  the  inverse  functions 

x = x(u,v)  y = y(u,v)  (2.34) 

are  single  valued.  For  multiple  values,  we  proceed  as  in  the  case  of 
a single  variable  and  obtain  a result  similar  to  (2.31).  The  last 
member  of  Equation  2.33  is  called  the  Jacobian  of  the  transformation 


3x 

3u 

|X 

3u 


3x 

3v 

3v 

3v 


(2.35) 


which  is  part  and  parcel  of  the  transformation.  We  must  caution  that 
failure  to  include  the  Jacobian  would  lead  to  erroneous  conclusions  re- 
garding PDF  of  the  new  variables. 

2.3  RANDOM  SIGNAL  WITH  UNIFORM  PHASE 

In  order  to  provide  a better  feel  for  the  material  or  probability 
theory  introduced  in  this  chapter,  let  us  consider,  for  example,  the 
ensemble  defined  by  the  sinusoids 


x(t,$)  - A cos(wt+(J>) 


(2.36) 


where  A,  to  are  fixed  and  $ distributed  uniformly  over  a complete  period. 
This  means  l 


P(4>) 


H* 


(2.37) 


In  practice,  x(t,$)  may  typically  represent  a scattered  field,  with 
A,  to,  and  $ representing  the  amplitude,  carrier  frequency,  and  phase, 
respectively.  The  various  probabilities  and  expected  values  can  be 
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determined  in  the  following  manner.  First  we  invert  Equation  2.36  to 
obtain  the  two  branches; 

-1/x 


4>x  = COS'1^)  -Wt,  (fK^fTT 
$2  = cos"*(^)  -Ut+TT,  TT<_4)2£2ir 


Di f f er entiat ing , 


d<j>2 

33 r 

33T 

(A2-*2) 


177 


so  that 


p(x)  = 


1 


.( 


A" -x“  ) 


,2  2 


T7T 


(2.38) 


Note  that  the  result  does  not  depend  on  to  or  t (stationary) , As 
x •*■  ±At  p(x)  becomes  infinite,  which  seems  co  conflict  with  the  funda- 
mental fact  that  probability  can  never  exceed  unity.  The  answer  lies 
in  realizing  that  p(x)  by  itself  has  no  physical  meaning  unless  it  is 
multiplied  by  the  differential  element  dx,  which  together  give  the 
probability  of  finding  x in  the  range  (x,x  + dx) . This  will  never 
exceed  unity.  He  conclude,  therefore,  that  the  signal  level  is  most 
likely  to  be  found  near  +A  or  -A,  If  A were  not  a constant,  one  has 


p(x?  , i_  | ■ --LW.dA 


(A2-2) 


m 


(2.39) 


where  p(A)  is  the  PDF  of  A. 


For  the  second  order  PDF,  all  that  is  needed  is  the  conditional 
probability  occurring  in  Equation  2.3.  Since  Equation  2.36  is  a 
deterministic  function,  once  its  value  is  known  at  t^,  it  is  speci- 
fied (functionally)  for  all  other  times.  Thus 
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X2  = A cos  ^ ujt 2+4> ) 


A cos 


OJT  + COS 


■‘(r) 


(2.40) 


with  T = as  usv..-, ' , 

Employing  the  delta  function  notation  to  denote  the  PDF  of  a 
constant/  we  have 


e(x2|V*)  - 6 


X2  - A COS  { U>T  + cos 


(“T  * "’’Hr)) 


(2.41) 


It  is  much  easier  to  compute  ensemble  and  time  averages,  for  example. 


E(x)  = / xp(<i>)d<$)  = 0 

-u 


(2.42) 


x * A cos  (u>t+<j>)  » 0 


(2.43) 


xlx2  3 T~  cos 


(2.44) 


E ) " r~  cos 


(2.45) 


We  note  in  particular,  the  process  is  not  only  stationary,  but  ergodic 
in  the  sense 


E(x)  ■ x 


E(xix2)  ■ Xlx2 


(2.46) 

(2.47) 
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If  the  restriction  that  A be  a constant  is  removed,  the  second 
relation  will  not  be  true  in  general,  which  means  that  the  process  fails 
to  be  ergodic  with  respect  to  the  auto-correlation  function.  Al&e  if  4> 
is  not  uniformly  distributed  the  process  is  no  longer  stationary.  In 
what  follows,  we  will  assume  a uniform  distribution  for  <p. 
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SECTION  III 

SCATTERING  FROM  CHAFF  CLOUDS 


3.1  GENERAL  CONSIDERATIONS  j 

The  word  chaff  denotes  a confusion  type  electronic  countermeasure  j 

employing  a large  number  of  resonant  dipoles.  The  dipoles  are  usually  1 

in  the  form  of  very  narrow  aluminum  strips  or  aluminum-coated  glass 
fibers  cut  to  a length  of  about  one-half  wavelength  at  the  frequency 
of  interest.  Since  the  bandwidth  of  narrow  dipoles  is  quite  small, 

one  generally  uses  several  cuts  to  obtain  coverage  over  a wide  band  jj 

of  frequencies.  When  properly  distributed  in  space,  a chaff  cloud  may  i 

l 

occupy  a large  volume.  At  microwave  frequencies  a single  chaff  pack- 
age contains  literally  hundreds  of  thousand  dipoles. 

{ 

I 

For  a distance  R from  the  transmitter  that  is  large  compared  to  j 

the  pulse  width  T,  the  number  of  scatterers  or  dipoles  per  resolution  I 

cell  is  ; 

I 

f 

N " n (Ifi)  2 |I  ,-3.!, 


where  n is  the  average  number  of  scatterers  per  unit  volume,  0 the 
radar  beamwidth,  and  c the  velocity  of  light.  At  any  instant  of 
time,  one  may  assume  that  the  number  of  scatterers  entering  the  rang? 
cell  is  equal  to  that  leaving  so  tnat  N can  be  considered  more  or  less 
constant  and  not  a random  quantity.  Making  N random  will  complicate 
the  problem  unnecessarily  without  altering  the  conclusions  in  any 
substantive  manner. 

It  is  obvious  that  a very  large  number  of  elementary  targets  are 
involved  in  the  scattering  process.  The  signal  scattered  by  each 
elementary  target  will  have  random  phase  and  amplitude  because: 

(1)  the  orientation  of  the  dipole  may  change  due  to  rotation,  and 

(2)  the  distance  between  radar  and  the  dipole  center  may  change. 

The  phase  of  the  returned  signal  is  actually  independent  of  the 
orientation,  but  the  amplitude  is  not.  However,  if  the  rotation 
rates  are  small  compared  to  time  of  observation,  the  amplitude 


changes  may  be  neglected.  The  dependence  of  return  power  on  Initial 
orientations  can  be  included  by  suitably  defined  PDF  for  the  angles. 

The  second  cause  can  affect  amplitude  as  well  as  phase.  The  amplitude 
changes  are  quite  small  and  may  be  ignored.  The  phase  changes  are 
most  important  and  have  been  studied  quite  thoroughly.  One  generally 
assumes  that  the  phase  distributions  are  uniform  over  a full  cycle. 

What  this  means  is  that  a single  dipole  may  occupy  any  position  within 
the  range  cell  with  equiprobability.  Also,  the  number  of  dipoles  wich 
any  given  phase  will  be  the  same  as  those  with  any  other  phase.  The 
dipoles  will  also  be  assumed  to  be  independently  moving  and  the  effects 
of  mutual  coupling  will  be  neglected.  The  effects  of  wind  will  be 
examined  to  some  extent.  Since  we  are  dealing  with  noncoherent 
scattering,  mass  motion  of  the  cloud  with  constant  speed  will  have 
negligible  effects  because  all  the  individual  dipoles  are  affected 
equally.  However,  if  there  is  relative  motion  between  the  dipoles 
either  due  to  turbulence  conditions  or  some  other  reason,  the  return 
signal  will  fluctuate  proportionately.  These  fluctuation  rates  are 
attributable  to  the  doppler  beats  of  the  individual  scatterers.  Using 
probabilistic  methods  we  will  now  develop  the  first  and  second  order 
statistics.  The  material  has  been  gathered  from  several  sources  shown 
under  references.  Some  of  the  derivations,  especially  second  order 
statistics,  are  obtained  by  simpler  means. 

3.2  FIRST  ORDER  STATISTICS 

For  a collection  of  N scatterers,  the  resultant  complex  signal  S 
is  given  by  the  vector  sum  of  the  individual  returns  and  if  one 
neglects  multiple  scattering, 

S - Ve10  - ? A(k)  ei4,Ck)  (3.2) 

k-1 

where  A ' is  the  amplitude  and  <J>  ' the  phase  of  the  k scatterer  and 

for  convenience  the  additional  phase  term  due  to  the  carrier  frequency 
has  been  factored  out.  Our  problem  is  then  to  find  the  probabilities 
of  V and  0 given  the  probabilities  of  A^k*  and  . The  above  sum 
represents  the  familiar  random  walk  problem  in  the  complex  plane. 
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Resolving  S into  real  and  imaginary  parts,  we  have 


w 


Re  S - V cos  8 - x - T A(k)  cos  4»(k)  (3.3a) 

k=l 


N 


.00  c -i  n *00 


Im  S = V sin  0 = y = [ A • * sin  <p 

k=l 


(3.3b) 


(k'  (k) 

Knowing  the  joint  PDF  of  A ' and  $ , one  might  compute  the  proba- 

bilities of  the  individual  terms  and  the  sums  by  the  methods  outlined 
in  Chapter  II.  However,  in  view  of  a very  powerful  theorem  called  the 
Central  Limit  Theorem,  thera  is  no  need  to  go  about  this  the  hard  wc;y. 


Central  Limit  Theorem: 

Let  X^,  X(2),...  be  N independent  random  variables  all  of  which  have 

2 

the  same  distribution  with  expectation  U and  variance  O . The  dis- 
tribution of  the  sum 

S = f X{U  (3.4) 

k*l 

2 

approximates  normal  for  large  N with  E(S)  * Ny  and  D(S)  » so  . 

In  other  words 


P(S) 


-(S-Nu)2 

J e”  2o‘N 

/TiF  a 


(3.5) 


1 i 

S 

\ 


The  beauty  of  the  theorem  is  that  one  need  not  know  or  be  concerned 

with  the  individual  distributions.  It  is  not  known  precisely  how  large 

N should  be,  but  the  conditions  are  almost  always  satisfied  for  chaff 

clouds.  Let  us  for  a moment  assume  that  all  t.ie  dipoles  are  cut  to 

the  same  length  (later  we  shall  account  for  multiple  length)  and 
(k)  (k) 

A , ♦ are  uncorrelated  which  is  justified  in  view  of  earlier  dis- 
cussions, as  a matter  of  fact  they  arc  statistically  independent.  Let 
the  distribution  for  phases  be 
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I 


P(4>) 


27 


, - 7T  < 4>  < 7T 


(3.6) 


Now 


eTa'^a'^ 


cos  <i>^cos  <f>^ 


T 


&i 

1 


6|  is  Kronecker  delta 


(3.7) 


and 


E[Vk>cos  ^k>]  - E [ACkJ]  E [cos  *(k>] 


= 0 

Invoking  the  central  limit  theorem 


- r 

P(X)  = — o 2o' 
/It r o 


(3.8) 


2 N 2 

where,  since  the  amplitudes  have  identical  PDF,  a = ~ E(A  ) and  a 
similar  expression  for  Y,  of  course.  To  proceed  further,  we  need  the 
joint  PDF  of  X and  Y.  In  general,  it  is  not  possible  to  write  down 
the  joint  density  from  a knowledge  of  marginal  densities.  An  excep- 
tion to  this  is  the  normal  process  where  all  order  densities  are 

normal.  For  the  case  of  two  variables,  with  zero  mean  and  identical 

2 

variances  o , the  result  is  particularly  simple, 


|HX,V)  a 


Iso" ( 1 -o' ) 


TH’  ex(  > 


r x---.3\v«v: 

L (i •?•) 


where  the  correlation  coefficient  is 

„ . MML. 

[ l:  ( X" ) t ( v” ) ] 


171 


(5.9) 
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Prom  Equation  3.3,  we  note  that  E(XY)  = 0,  so  that 


2 2 
-X  +Y^ 


2a 


P(X,Y)  = 


2na‘ 


« p(X)p(Y) 


Introducing  polar  coordinates 


1/2 


V = (X2+Y2)  , 8 = tan'1  ( £) 


we  have,  using  Equation  2.33 


p(V,0)  = 


V 0 ^ 


and 


IQ  = 2o2  = NE  (a2  ) 


(3.10) 


(3.11) 


The  marginal  densities  are  obtained  by  integration  as  usual. 


P(V) 


V2 

S.’  *» 

‘o 

o < V < » 


(3.12) 


1 


P(6)  “ J7  • •H<(M  8 


(3.13) 


„I 


. f ' 
•il  .' 


The  resultant  phase  is  uniform  rad  the  distribution  for  V is  called 
Rayleigh  or  x with  two  degrees  of  freedom.  These  results  have  been 
known  for  about  BO  years  and  apply  to  many  kinds  of  scattering  situa- 
tions. The  amplitude  distribution  of  the  elementary  signals  does  not 


enter  into  the  picture  and  it  is  therefore  quite  erroneous  to  assume 
that  the  dipoles  in  a cloud  should  have  spherically  uniform  orienta- 
tions to  arrive  at  the  above  statistics.  The  requirement  of  uniform 
phase  is  essential,  however.  For  nonuniform  phase  distributions,  the 
analysis  is  very  similar,  but  the  results  become  quite  complicated. 
The  important  expectations  in  the  present  case  are 


/it  I 

E(V)  = -y-2- 

E (v2)  ' i0 

E(0)  - 0 
E (e2)  « tt2/3 


(3.14) 


(3.13) 


The  cumulative  distribution  of  V is 


•v2/i 

P(V)  - 1-e  0 


(3.16) 


From  a practical  viewpoint  the  radar  cross  section  (RCs)  of  the 
2 

cloud,  defined  by  I * V is  more  significant.  One  has 

’UK 


p<») 


, I **0 


(3.17) 


and 


i’(U  « l e 


1/1. 


is.  m 


2 


The  distribution  is  called  Rayleigh  power  or  x with  two  degrees  of 
freedom.  It  is  easily  checked  that 
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E(I) 

E (i2) 


= I. 


= 21: 


SD  = I 


(3.19) 


One  normally  expresses  the  SD  as  a percentage  of  the  expectation  which 
in  the  present  case  is  100%.  A plot  of  the  distribution  of  V and  I are 
shown  in  Figures  1 and  2.  (Observe  that  the  most  probable  value  of  RCS 
is  zero.)  The  median  is  0.7  IQ  and  the  signal  intensity  will  be  one 
half  the  average  level  39%  of  the  time. 


The  results  obtained  thus  far  are  quite  general  and  let  us  now 
specialize  to  chaff  clouds.  All  that  remains  to  be  done  is  to  relate 
Iq  to  the  dipole  scattering  properties.  For  a half  wave  dipole  in- 
clined at  angles  0 and  <p  in  a spherical  coordinate  system,  the  back 
scattered  amplitude  A for  a plane  wave  traveling  in  the  z direction 
can  be  approximated  by 


2 2 

A * /a~  sin  0 cos  <j> 
o 


(3.20) 


where  in  terms  of  wavelength  A,  the  broadside  RCS  o = 0.89A  . 

o 


What  we  now  need  is  the  joint  PDF  p(0,<J>)  describing  the  orienta- 
tions of  the  dipoles.  It  depends  upon  the  rotation  rates,  emironment, 
etc.,  in  a complex  fashion  and  no  serious  attempt  has  ever  been  made 
in  this  direction.  One  therefore  makes  the  simplifying  a priori 
assumption  that  all  orientations  are  equally  likely  so  that 


p(e,)>) 


sin  0 0<_<)k27t 


4 TT 


(3.21) 


0 < 0 < IT 


whence 


IQ  = NE  (a2)  = 


No. 


(3.22) 
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This  is  the  so-called  tumble  average  RCS  of  the  chaff  cloud.  In  case 
of  orthogonal  polarization  reception,  the  above  expression  should 
be  divided  by  three.  For  chaff  clouds  that  contain  several  cuts  of 
dipoles,  we  simply  take 

1Q  - N1o1  + N 2a2  +...  (3.23) 

where  the  N's  refer  to  number  and  0's  suitably  defined  averages. 

3.2.1  Steady  Target  Immersed  In  Chaff 

The  determination  of  the  radar  return  of  an  aircraft  flying 

through  a chaff  corridor  is  of  considerable  practical  importance.  The 

radar  cross  section  of  an  aircraft  changes  with  aspect  considerably  and 

there  are  no  reliable  statistics  on  the  subject.  We  will  therefore 

consider  the  aircraft  as  a steady  target  and  determine  the  combined 

statistics  of  aircraft  plus  chaff.  This  will  be  at  least  first  step 

toward  understanding  a more  difficult  problem.  If  the  RCS  of  the 
2 

steady  target  is  m I , we  have  for  the  probability  densities 


P(V,0) 


Ve 


/ 2 V ~ v -2Vm 
- ( m + v—  1 cos  0 

V lo  ) /T 

e 0 


(3.24) 


Integrating  over  0 


P(V) 


2Ve 


'("2*0 
T 


(■?) 


(3.25) 


and  also 


Pd) 


• •'  w (••■  /t) 


(3.26) 
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where  J is  zeroth  order  Bessel  function.  The  distribution  for  phase 
o 

is  quite  complicated  and  will  not  be  given  here.  The  above  is  called 
the  Rice-Nakagami  distribution  originally  discovered  in  the  study  of 
fading  of  radio  signals  due  to  mixed  propagation  paths.  Using  the 
known  integrals, 


(3.27) 


,8/ot 


2a2  + 4a£  + B2 


(3.28) 


= f0  X2e  ax  Jq  (i2/Bx)dx 


we  get 


E(I)  ■ (l*»2)(l0) 


(3.29) 


which  could  have  been  guessed  and 


E (l2)  - I2  (2  + 4m2  ♦ m4) 


(3.30) 


SD  * I (l  + 2ra2) 


(3.31) 


i 


SD  . (l  ♦ 2m2) 

1 ♦".* 


(3.32) 


5 /I/m,  for  m2>>l 


if.' 


- *.‘-v  -.1 


Summarizing  we  note  that  RCS  statistics  of  a chaff  cloud  with 

,2 


uniform  distribution  of  dipoles  follow  X distribution  with  two  degrees 
of  freedom.  If  there  is  a nonfluctuating  target  in  addition,  we  have  a 
Rice-Nakagami  distribution.  The  fluctuations  are  quite  large  (100%) 
in  the  first  case  and  the  presence  of  a strong  steady  target  decreases 
the  fluctuations. 


3.3  SECOND  ORDER  STATISTICS 

Whereas  the  first  order  statistics  tell  us  about  the  magnitude 
of  fluctuations,  information  regarding  the  rates  of  fluctuations  is 
obtained  from  the  second  order  statistics  which  give  the  probability 
of  jointly  finding  two  values  of  a random  variable  at  different  times. 
Analogous  to  Equation  3.3,  we  define  the  joint  density 


p ^X1,X2;t1,Y1>Y2;t2^ 


(3.33) 


where  X^  * X(t1> , X2  = X2(t2),  etc.  First,  we  assume  the  process  is 
stationary  so  that  time  appears  as  only  the  difference  T * t^t^. 
Also  since  X and  Y are  statistically  independent,  we  have 


P (xi*X2;Y1,Y2;t)  - p (X^X^t)  p (Y^Y^t)  (3.34) 


We  first  observe  that  the  marginal  densities  of  X'.  and  X.  are 


normal  and  thus  given  by  Equation  3.8.  Therefore,  analogous  to  Equation 
3.9  we  have 


exp 


p (Xj,X2; 


f 1 

L I„  (1-62)  j 

I7T 


»i0  (i-s  ) 


(3.35) 
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to  obtain 


V1V2 


p (Vj.Vj.ej.e^t)  = 


exp  - 


rv^v|-2*v1v2  cos(ere2) 

■„  (l-s2) 


o<V^,V 2<0&»  “7rl9^>02iir 


Since  the  above  cannot  be  written  in  the  product  form 


p ( V1 * V2 ) P (01 *62  ) 


the  amplitude  v and  phase  6 are  correlated. 


Integrating  over  the  angles 


p (V1,V2't) 
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L' TTT^T  J 


exp 
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Utt^J 


•As 


(3.39) 
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If  we  denote  the  ratio  V^ | V^  by  Q 


p(q,v2;t)  - 


4QV: 


T(w2) '° 


2gQV2  1 

1 ^T^g1)  J 


exp  - 


r v2 
v2 


( 1+Q2  ) 

X 


>0  (*V) 


and 


P(Qit)  = / p (Q,V2;t)  dv2 


using  the  integral 


/ Xe  aX  J (iBX)dx  ° 
o 


(a2-*2) 


T7T 


the  distribution  is  found  to  be 


p(Q;t) 


20  1 

(.•g2)  1 

(l+Q2 

) 

[ (l*Q2) 

i -w] 
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(3.41) 


(3.42) 


(3.43) 


(3.44) 
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s'  i'v; 
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and  is  plotted  in  Figure  3.  Note  that  as  g -*■  1,  it  peaks  at  the  center 
like  a delta  function.  Similarly  one  has  for  the  intensities 


and  for  the  ratio  I2^Ii'  denoted  w» 


p(W;x) 


('-E2) 

7 22 
[(1+W)  J 


0<w<‘ 


(3.46) 


This  distribution  is  shown  in  Figure  4.  To  get  the  joint  density  for 
6,  we  integrate  Equation  3.39  over  and  V2.  First  we  introduce 
polar-coordinates 

Vj  - p cos  <f>,  V2  ■ p sin 


so  that 


r ~ ' 'Pf Tif  rumiif i n i -•-.:  "■-"  ■:.  * .'-. '.  2iL Si  -"'; 


•yy i r'Wf^Cfiw 


sssIwSm^t*??^™?3^^  ■ ‘:'ryy  • 

1 ' V.  ~ " ' ihm unii'iiifrtri  'SM^  y 


P (9l»e2;T) 


tt/2 


172  77  2 v / sin  2*  d* 


2^1^  (l-g2) 


/ p-5  exp 

o 


2 1 
p (1 -y  sin  2^) 


'.M 


dp 


with 


y 3 g cos  ( 9i *®2  ) 


Noting  that 


/"  XV“2,(2  dx  - 1 
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the  above  becooes 


1 -|2  n j 2 sin  2i|» 


4n  0 ( 1 - y sin  2t|») 


di|< 


1- 


.2  it 


sin  X 


8«‘ 


(1 *y  sin  X) 2 


dX 


1-K2  d J*  dX 


8k  * 


3y  q (1 • y sin  X) 


?9 


»• 


\ 


'»m*rn!SK 


and  using  the  result 


w a finally  have 


o • -1 

it  + 2 sin  y 


(>V) 


P (V02:t)  - 


(»V) 


ki!  fj_ 

4it  1-' 


/ n -1  ^ 

\2  + Sin  ^ 


(3.47) 


where 


Y B g(t)  COS  (Vh) 


-ti<_9  j *82!^ 


A more  convenient  expression  is  obtained  by  introducing  a new  set 


of  variables 


u ■ 0 ^ - 0 2 , v ■ Gj»  with  the  limits 


-»<u*v<3  an'i  -3<v<» 


The  total  probability  for  u is 


P(u;t)  = /*  p(u,v;-r)dv 


where  care  must  be  exercised  in  choosing  the  limits  of  integration 
Manipulating  the  inequalities  for  u and  v,  we  note 


a = Max(-ir, -it-u) 
8 = Min(it,ir-u) 


which  leads  to  two  possibilities 


u>o,  a = -h,  8 = *-u 
u<o,  a = -Ti-u*  8 = " 


Hence 


P(u;t)  - i’4-  (2t?  - | 
4ti* 


K cos  u 


U|)  r - — * — 2 — - 

(l-g*  cos  * u) 


( 1 -g2  cos2  u ) 


'mil*  sin  (8  cos 
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t?'> ; 
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f' 


IP 
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A very  useful  alternate  form  due  to  Middleton  is 


P(u;t)  = ^4hJ 
4 u 


1 ♦ 8 35  F(u) 


where 


F(u) 


+ sin  *(g  cos  u)]  sin  u 
= m 


(l-g: 


cos 


u) 


The  PDF  of  u is  plotted  in  Figure  5.  This  completes  the  derivation 
of  first  and  second  order  statistics  for  the  amplitude,  intensity 
(RCS)  and  phase  of  the  scattered  signal  by  a chaff  cloud. 
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SECTION  IV 

CORRELATION  FUNCTIONS 


4.1  CLOUD  CORRELATION  FUNCTION  g(r} 

The  correlation  function  g(x)  defined  by  Equation  3.36  is  about  the 
most  significant  quantity  characterizing  the  second  order  statistics. 

It  is  governed  by  the  dynamics  of  the  chaff  cloud  itself  and  only  in 
special  cases  is  it  possible  to  establish  the  connection  between  the 
movement  of  chaff  dipoles  and  g(i).  Now  if  A^k) , ^k)  and  A(k),  <j>(k) 
denote  the  amplitude  and  phase  of  the  scatterer  at  times  t and 
t2  ~ fcl  + T'  then  one  1133  for  the  total  field 


X,  = 


N 

I 

k=l 


A<«  cos 


(4.1) 


N 

I 

k=l 


cos 


(4.2) 


Simple  calculations  show  that 


E(X12)  "IE(A!2) 

H (XX)  > § B(a\) 


assuming  identical  scatterers.  Furthermore 


N N 


(X1X2  ) “ £ E (All)A2j))  E (cos  *x(l)  cos  (4.5) 


Since  the  scatterers  are  independent 

B (cos  cos  ) ■ 0,  i j j 
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contributions  occurs  only  for  i = j.  Therefore 


g(t) 


2E  (^^2)  E ^cos  ^ cos  ^2) 

[E(A0E(A2)]  1/2 


= 2E  ^cos  cos  $2  j 

(4.6) 


if  we  assume  2^  = which  is  in  keeping  with  the  understanding  that 
amplitude  changes  are  not  significant.  In  terms  of  the  change  in 
phase  between  the  two  instants  of  time,  let 


so  that 


l|)  ~ $2  ■ <j>l  , -“<t|)<00 


g(x)  = 2E  cos  <S>1  cos  (<i'+4>1)  j 


E(cos  tji) 


(4.7) 


(4.8) 


kn  identical  result  will  be  obtained  by  considering  the  Y-component. 
If  a dipole  moved  a distance  £ in  the  direction  of  the  radar  during 
the  time  interval  x,  then 


if* 


4 TI 


(4.9) 


where  X is  the  wavelength.  If  all  directions  of  motion  are  equally 
likely  arising  from,  say,  an  isotropic  turbulent  wind  field,  then 
the  joint  density  is 


P(V,0,$)  - sin  Oq(V) 


0<V<®,  0£8£n , o<$<2« 


(4.10) 
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where  V is  the  speed  and  8,  $ suitably  defined  angles  with  0*0  de- 
noting direction  away  from  the  radar.  Also  note  that 


/ q(V)dV  = 1 
o 


(4.11) 


Since  £ = VT  cos  8,  we  have 


PU,V,8,$)  = ptt|v,e,*)  POM,*) 


= 6(C-Vx  cos  6)p(V,e,$)  (4.12) 


and  therefore 


d£  n 00 

PU)<U  * y~  j sin  tidQ  J q(V)6(S-Vt  cos  0)dV  (4.13) 


Following  Sieyert, ^ we  express  the  delta  function  as  a Fourier 
integral 


6(S*Vt  cos  8)  » ^ J09  d ueiu^*VT  Cos  0)  (4.14) 


so  that,  after  integrating  over  8, 


P(C)  ■ t / dueluC  / dVq(V)  (4.15) 

-*  o 
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and  using  Equation  4.9 


pO) 


_ A 


8n 


00  W-  CO 

/ dpe^  / dVq(V)  iiSMll 
-00  o Mv  1 


and 


00 

g(t)  = E(cos  i JO  = / cos  ^ p(<|))dij) 

-00 


The  integrations  over  u and  tp  can  be  carried  out  in  the  following 
manner. 


iuXijj 


“2-  / / dydij>  cos  e 471  ^7'^Vt} 


7 //  cos  * ^^d*dU 


7? 


J*dusi^Vll  /*„,  [cos(l  .jji)* 

* COS  I 1 - ^4  \ l|»  1 


and  finally 


g(t)  - / q(V) 
o 


sin(4irVT/X] 

C4  ivxhy 


(4.1£ 


We  may  now  make  the  following  observations.  If  a chaff  cloud  is 
measured  simultaneously  at  different  frequencies  and  g(T)  is  plotted 
against  X/y  all  the  curves  should  coincide.  From  a frequency  domain 
viewpoint,  it  means  the  doppler  beats  are  proportional  to  the  carrier 
frequency  which  is  a well  known  fundamental  result.  By  inverting 
Equation  4.19 


2 » 


q(V)»|^(~)  / Tg(i)sin(4TrVx/X)dT  (4.2C 


which  may  be  utilized  to  estimate  the  speed  distribution  from  the  mea 
sured  value  of  g(t). 

4.2  SIGNAL  CORRELATION  FUNCTIONS 


Using  the  definition  (Equation  2.16),  we  have  for  tho  intensity 


(RCS) 


B(I,t)  - //  lxl2  pCIj.I^tJdljdlj 


("*2) 


(4.21 


and  the  auto-covariance 


This  is  an  interesting  result  because  it  provides  an  indirect  method 
of  determination  g(T)  from  the  RCS  auto-correlation  function.  The 
latter  can  be  determined  experimentally  using  a noncoherent  pulse 
radar.  For  the  amplitude  we  have 


B (V , x ) = //  VXV2  P(V1,V2;x)dV1dV2 


using  Equation  3.40,  and  introducing  new  variables 


Vx  - t1/2  cos(<j>/2),  V2  = t1/2  sin(4>/2) 


we  have 


2 2 1/2 

by  use  of  Equation  2.43,  and  letting  A - (1  - g sin  $)  . 

Using  elliptic  integrals 

I r _ 

B(V,t)  - [ 2^U)  * (l*g2)K(g)J  (4.23) 


\ 
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where  K and  E are  complete  elliptic  integrals  of  the  first  and  second 
kind  with  modulus  g.  Furthermore,  the  auto-covariance  is 


R(V,t) 


= 2Ete)  - 


(l-82) 


2 -it/2 


mi 


To  determine  B(8,t)  consider  the  second  moment  of  u = 


E(a2)  - E ^(6,-63)  2J 

• 2E  ( e2)  -2B(e1e2) 

* ~j*“  * 2B(0 ,t) 

from  which 

b(6,t)  - ji  - 4u2) 

Now  from  Equation  3.49 

E(u2)  » / u2p(u;t)c’.u 

' 7 -2* 

■ r /’  ox(x2-  *1)  [ i*8  FCx)  ] 
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(4.24) 


(4.25) 


(4.26) 


(4.27) 


The  final  integrations  can  be  carried  out  by  expanding  the  inverse 
sine  functions  into  power  series  and  integrating  term  by  term.  First  we 
note  below  the  known  results 


sin  (g  cos  X)  = 


- Z 


(2k) 1 (g  cos  X) 
~ ■ 2 k , | \ 2 


2k+l 
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Then,  it  follows  that 


2ir  _ , 

J sin'  (g  cos  X)dX 


2tt  _ -j 

f X sin  (g  cos  X)dX  = 0 
o 


/ sin'1  (g  cos  X)  2dX  = 7T  T]  I 
o n 

n=l  n 


2n 


2w  2 

/ X sin’* (g  cos  X)  2dx  * ^ ^ 


g2" 

iT- 


n=l 


Substituting  the  above  result  in  Equation  4.29  we  have  finally 


2tr 
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. -1 


’ (sin'1  g) 


E ( u ) * -y-  - tr  sin  g \sin  * g, 

nr> 

2n 


4 £ 


n°  1 


(4 . 30  j 


and 


E (V1’)  “ I [ 


it  sin'1  g ♦ 


(sin'1  g ) 


! * 2n  "I 

lj  z V 

n- 1 n J 


(4.31) 
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The  normalized  phase  auto-covariance  becomes 


R(0,t) 


00 


n=l 


(4.32) 


4.3  NUMERICAL  RESULTS 

The  problem  of  determining  the  speed  distribution  of  the  dipoles 
in  a chaff  cloud  under  a given  set  of  environmental  conditions  is 
extremely  complex  and  so  far  no  serious  attempts  have  ever  been  made 
in  this  direction.  We  will  therefore  be  content  by  presenting  results 
for  an  assumed  distribution.  For  example,  if  all  the  dipoles  have 
the  same  speed  Vq, 

q(V)  = 6 (V-V0)  (4.33) 


and  from  Equation  4.19 


sin  (4nV  t/A) 

«“>  ■ Ayr-- ■ 

The  auto-covariance  function  for  the  intensity  becomes 


(4.34) 


R(I 


sin  ( 4 w V t/a}  1 " 
t%)  & — 


and  the  frequency  spectrum  is 
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(4.35) 


(4.3b) 


0 


The  spectrum  thus  contains  frequencies  up  to  4 Vq/^  which  is  no 
more  than  the  doppler  beat  between  dipoles  moving  directly  into  and  away 
from  the  radar  beam.  It  should  be  emphasized  that  we  are  talking  of 
frequency  spread  due  to,  for  example,  turbulence  and  not  the  conven- 
tional doppler  frequency  due  to  average  motion  of  the  entire  cloud. 

The  latter  cannot  be  measured  by  a noncoherent  system  because  the 
phase  information  is  lost.  For  the  special  value  of  g(T)  given  by 
Equation  4.34,  the  three  auto-covariance  functions  are  plotted  in 
Figure  6 and  Figure  7 shows  typical  experimental  results  obtained  in 
a recent  AFAL  contractual  effort. ^ Although  the  data  in  the  two 
figures  are  unrelated,  one  notices  certain  trends.  For  the  type  of 
chaff  payloads  employed,  the  correlation  times  were  found  to  be  in  the 
order  of  10-20  milliseconds  which  means  that  the  frequency  components 
are  around  50-100  Hz. 


45 


I 


w 

ft 

U'-: 


I* 


m 


m. 


m 


m 


i 


i 


i 


m 


■r- 

I: 


% 


■ FREQ.  BAND. 
■TRANS.  POL. 


■ REC.  POL. 


0.2  03  0.4 

LAG  TIME  X FREQUENCY  {QHi-Sk.I 


FREQ.  BANDS 


E - 2.8GHz 
G - 5.5 
I - 9.25 


POL.  SYMBOLS 


H - HORIZ. 

V - VERT. 

R - RIGHT CIR. 
L - LEFT  CIR. 


Figure  7 


Intensity  Auto-correlation  Function  by 
Experiment 
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SECTION  V 

CONCLUSIONS  AND  RECOMMENDATIONS 


This  effort  accomplished  the  task  of  obtaining  mathematical 
expressions  for  the  first  and  second  order  statistics  of  electro- 
magnetic scattering  from  chaff  clouds.  These  results  can  form  a basis 
for  studying  the  effects  of  chaff  clutter  fluctuations  on  different 
types  of  advanced  radars  in  a statistical  sense.  For  example,  the 
starting  point  in  the  case  of  noncoherent  MTI  with  a single  delay  line 
canceller  is  the  conditional  probability  for  the  amplitudes  given  by 


P(V2|Vt) 


P ( V1 ,V2 ’ T ) 


P /V 


(M 


2VZ 

= _ . T 

r . w2  i 

Xo  ( 

>•*)  ° 

J 

exp  - 

V2*B2V2 

Jo  (W). 

(S.ij 


where  i is  the  interpulse  period  and  and  are  the  amplitudes 
for  two  successive  pulses.  The  uncancelled  output  is  [V^-Vj. 

Herein  lies  the  physical  significance  of  q.  If  g *=  1,  there  will 
be  perfect  cancellation.  A quantitative  study  can  now  be  made  in  terms 
of  the  correlation  functions  of  chaff  and  the  target  as  parameters. 
Details  of  this  and  similar  studies  will  be  presented  in  subsequent 
reports. 
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